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Abstract 

Currents in a metallic ring with a quantum dot are studied in the framework of a Langevin equation for a magnetic flux passing 
through the ring. Two scenarios are considered: one in which thermal fluctuations of the dissipative part of the current are modelled 
by classical Johnson-Nyquist noise and one in which quantum character of thermal fluctuations is taken into account in terms of a 
quantum Smoluchowski equation. The impact of the amplitude and phase of the transmission coefficient of the electron through 
a quantum dot on current characteristics is analyzed. In tailored parameter regimes, both scenarios can exhibit the transition from 
para- to diamagnetic response of the ring current versus external magnetic flux. 
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1. Introduction 

In the early 90’s after the successful reduction of the signal- 
to-noise ratio the three groups conducted pioneering experi¬ 
ments with the mesoscopic metallic rings. The careful mea¬ 
surements of Cooper (T), Gold (2), and Gallium-Aluminum- 
Arsenide/Gallium-Arsenide 0 normal rings has shown the ev¬ 
idence of the existence of the persistent equilibrium currents 
flowing in the small metallic pieces of the rotational symmetry 
reaffirming an old idea of Friedrich Hund (4). This very idea 
concerns the charge transport in normal metallic ring. From 
the Ohm law we can expect that from the macroscopic point 
of view such current will die out within the relaxation time for 
a given material, which for metal is known to be rather short 
and of the order of lO -14 ^. However, for sufficiently small cir¬ 
cumferences the macroscopic description is no longer valid and 
ring reaches the region where both macro- and micro-world 
meet making the requirement for the mesoscopic description 
0 of the dynamics. In low enough temperature the effects of 
quantum coherence of electrons appear. Under the right cir¬ 
cumstances some electrons in the ring are able to preserve its 
coherence which in turn results in a persistent (dissipationless) 
equilibrium current induced by the static magnetic field. In 
1965 Bloch 0 and five years later Kulik 0 confirmed Hund’s 
theory using the quantum-mechanical description. The real in¬ 
terest in the topic of the persistent currents in normal rings arose 
after 1983 paper by Buttiker et. al [8] where the existence of the 
persistent currents was shown also in the presence of the elastic 
dispersions. 

First measurements of currents in the diffusive regime 0 
have shown rather strong disagreement (10-200 times larger 
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currents amplitudes) with the theoretically anticipated values. 
Later attempts reduced this dissimilarity to a factor of around 2- 
3 0. Experiments with semiconducting materials in the close 
to ballistic regime usually agreed with the theory mm. Only 
recently the scanning SQUID technique was used to record not 
only the response signal of the rings itself but also from the 
background. This method gave the possibility of the high pre¬ 
cision measurements of the current flowing in 33 different sep¬ 
arate Gold rings ini and finally confirm qualitatively as well 
as quantitatively all aspect of the existing theoretical descrip¬ 
tions tm The alternative method was used to measure the 
currents in the Aluminum rings which were deposited on a can¬ 
tilever da. A torque magnetometer whose vibration frequency 
can be precisely monitored was used as a detector. The mea¬ 
surements was performed with the several different cantilevers 
decorated with a single aluminum ring or arrays of hundreds 
or thousands of identical Aluminum rings. The analysis of the 
different magnetic susceptibilities seen in ITTT11T31 based on the 
two-fluid model was addressed in da da da. 



Figure 1: Schematic picture of the mesoscopic non-superconducting ring with 
a quantum dot QD. The ring is threaded by a magnetic flux 0. 
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the mesoscopic metallic (non-superconducting) ring with the 
quantum dot. The experiment with the measurements of the 
phase of the transmission coefficient through a quantum dot 
in the Coulomb regime was performed in 1995 fT71 . Many 
different aspect of the persistent currents in the same sce¬ 
nario was studied rather intensively over the last two decades 
Similar schemes with the mesoscopic ring 
coupled to the quantum dot (23] l24l (25] or the quantum ring 
surrounding the quantum dot - a dot-ring nanostructure (DRN) 
was also addressed. Here we follow the model 
proposed by Moskalets (29) for a mesoscopic ring containing a 
potential barrier with a resonant level. 

The work is organized in the following way: In Sec. 2, 
the model is described and the Langevin equation for the mag¬ 
netic flux is presented both in the classical and quantum Smolu- 
chowski regimes. Discussion of the results is presented in Sec. 
3. In Sec. the stationary probability distribution of the mag¬ 
netic flux is analyzed. In Sec. |3.2| the impact of parameters of 
the quantum dot on average stationary currents is studied and 
regimes of paramagnetic and diamagnetic response are worked 
out. Sec. 4 contains summary and conclusions. 

2. Flux dynamics of mesoscopic metalic rings with a quan¬ 
tum dot 

We consider a mesoscopic metallic ring in an external mag¬ 
netic field B e applied perpendicular to the plane of the ring. 
At zero temperature the ring can display a persistent current I P 
when the size of the ring is reduced to the scale of the elec¬ 
tron quantum phase coherence length and the thermal length. 
At non-zero temperature T > 0, a part of electrons loses phase 
coherence due to thermal fluctuations and this part of electrons 
contributes to a dissipative Ohmic current I R associated with 
the resistance R of the metallic ring. The total magnetic flux 0 
piercing the ring is a sum of the external flux (p e oc B e and the 
flux due to the flow of the current /, namely, 

0 — + LI. ( 1 ) 

Here, L stands for the self-inductance of the ring. The current I 
is a sum of the persistent and dissipative currents, 

I = Ip + I R . (2) 

Now, following Ref. (29), we assume that the ring contains 
a potential barrier with a resonant level (a quantum dot). The 
expression for the persistent current I P = 4(0) in such a system 
takes the form [29] 

oo 

Ip = Io G{<pl<po) ^ A n (T/T*) cos[n(k F l + 5^)] 

n= 1 

x sin4arccos[tp cos(27r0/0o)]}. (3) 

The flux quantum 0o = h/e is the ratio of the Planck constant 
h and the charge e of the electron, 7 0 is the maximal persistent 
current at zero temperature for the ring without the quantum dot 
and 

mix \ t F sm(2n<f>/<f>o) 

G(<p/<p o) = . (4) 

y 1 -t 2 F cos 2 (2n(f>/(f>o) 


modifies the maximal current due to the quantum dot. Here, t F 
and S F are the amplitude and phase of the transmission coeffi¬ 
cient Tk = 4 expfzYy through a quantum dot for an electron of 
the Fermi energy. The amplitudes A n (T/T*) are determined by 
the relation 


A n (T/r ) 


T/T* 

sinh (nT/T*) 9 


(5) 


where T* is the characteristic temperature which measures the 
level spacing at the Fermi surface. The magnetic flux 0 is quan¬ 
tized with the flux quantum 0o = h/e being the ratio of the 
Planck constant h and the electron charge e. Moreover, k F is 
the Fermi momentum and l is the circumference of the ring. 
For t F — 1 and 5/7 = 0 this expression reduces to a current for a 
pure metalic ring (30). 

According to Ohm’s law and Lenz’s rule, the dissipative cur¬ 
rent Ir - 4(0) assumes the form 


1 d(p 

R ~ ~R~dt + 



( 6 ) 


where k B is the Boltzmann constant. It means that we include 
the effect of a nonzero temperature T > 0 by adding Johnson- 
Nyquist noise T(t) which represents thermal fluctuations. They 
are modeled by ^-correlated Gaussian white noise of zero mean 
and unit intensity, 


<r(0> = o, <r(or(j)> = 6{t - s). (7) 

Inserting Eqs. ([6]) and ^ to the relation (1) yields 

i d ~dt = _ z (0 ■ ^ + ip(<p)+ V¥- (8> 

We note that this equation is a Langevin equation for the mag¬ 
netic flux 0 = 0(0. Indeed, it has the same form as a Langevin 
equation for an overdamped motion of a classical Brownian par¬ 
ticle subject to the force F = F(0) which reads 

F(0) = -i(0-0 e ) + 4(0) (9) 

and the noise intensity strength D = k B T/R is in accordance 
with the classical fluctuation-dissipation theorem El E2). 
Therefore we can apply the well-known mathematical and nu¬ 
merical methods for analysis of Eq. ([8]). First, we transform it 
to the dimensionless form (see (33][34) for details) 

£ = -^ + V2D^». (10) 

as ax 

where x = 0/0o is the dimensionless magnetic flux. The new 
time s = t/ro , where the characteristic time To = L/R. The 
thermal noise intensity D 0 = k B T/((p^/L) = {E t */E^Tq = 
koTo, where the dimensionless temperature To = T/T*, E p * = 
k B T*/2 is energy of thermal fluctuations at the characteristic 
temperature T*, E^ = 0q/ 2L is the elementary magnetic energy 
and ko = E T * /E# rescales intensity of thermal noise. Rescaled 
Gaussian white noise £(s) has exactly the same statistical prop¬ 
erties as the dimensional version T(t). The rescaled potential 
V(x) takes the form 

= (11) 


2 










The rescaled external magnetic flux is denoted by = 0 e /0 0 
and the nonlinearity parameter a = L/o/27r0o- The potential 
consists of the harmonic part (x - x e ) 2 1 2 and the periodic part 


WO c) 


y A n (T 0 ) 
n 


cos (nSp) 
x cos{ft arccos[^ cos(27rv)]}, 


and the problem in a general case is still unsolved. In the so 
called quantum Smoluchowski limit, the leading quantum cor¬ 
rections are incorporated in the modified diffusion coefficient 
Do (361 El M E3 |40l SH 02 S3 SI . The modified diffusion 
coefficient takes the form ED 


n= 1 


( 12 ) 


DJx) = 


1 


PH -AfiV"(x))’ 


P~ l = D 0 . 


(19) 


where Sp = kpl + dp is a shifted phase. 

The Fokker-Planck equation corresponding to the Langevin 


j, p{x • 

where P(x, t) is a probability density of the process determined 
by Eq. ©• From this equation, all statistical properties of 
the magnetic flux can be obtained. In particular, its statistical 
moments ( x k (t )) are determined by the expression 


d 

r dv(x) i 

d 2 


A = A 0 

7 + y(i + JL) 

dx 

-j—-P(x,t) 

dx 

+ D 0 —Pix,t), 

(13) 
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The prime denotes the differentiation with respect to x. The 
dimensionless quantum correction parameter 


( 20 ) 


where 


%R 

= ~T’ 
7T0 O 


n 
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2nCR k B T* 


( 21 ) 


X oo 

/ P(x, i)dx, k = 1,2,3,... 

oo 


(14) 


For experimentalists, more interesting is the electrical current 
flowing in the ring. From Eq. 0 it follows that at any time the 
total current reads 


I(t) = 1(0(0 - 4> e ) 

and its average value is given by the relation 

L 

lit) = (x(t)) - X e , lit) = —(lit)), 


The psi function 'F is the digamma function (the logarithmic 
derivative of the gamma function). The y « 0.577 is the Euler 
constant and C is capacitance of the system related to the charg¬ 
ing effects. The quantum correction parameter A is a difference 
between the quantum (x 2 ) q and classical ( x 2 ) c second statistical 


( 22 ) 


moments of the magnetic flux (see Eq. (5) in Ref. 1 136 1), 

A = (X 2 )a ~ (X 2 ) c - 


(15) 


(16) 


The modification of the diffusion coefficient © results in 
modification of the Langevin equation, namely, 


dx 


_ _dVix) + j2D A (x) g(s) 
ds dx 


(23) 


where the dimensionless current i{t) has been introduced. In the 
stationary state, 


i = (*) - x e , (x) : 


S^OO 


d 

d 

\dVix) 1 

1 v P(x)dx , 

J —oo 

(17) 

^(*, t) = 
dt 

dx 

P(x, t) 
dx 


and should be interpreted in the Ito sense E3. The correspond¬ 
ing Fokker-Planck equation has the form 


+ — [D A (x)P(x,t)]. (24) 


where P(x) = lim^oo P(x, t) is a stationary probability density. 
It can easily be calculated from Eq. © for dP(x , t)/dt = 0 and 
zero stationary probability current yielding the distribution 

P(x) = lim P(x, t) = No exp [-T'cMl (18) 

t^oo 

and No is the normalization constant. The generalized thermo¬ 
dynamic potential T'cM = V(x)/Dq depends on the external 
flux and the stationary probability density is given by the 
Boltzmann distribution. Eqs. ([T7])-([T8]) form a closed set from 
which the non-linear function i = f(x e ) can be calculated de¬ 
termining the stationary current-flux characteristics. 

2.7. Quantum Smoluchowski limit 

Thermal fluctuations modeled as classical d-correlated white 
noise are adequate to describe many physical phenomena even 
in low temperatures. However, in some low temperature 
regimes, quantum effects like tunnelling, quantum reflections 
and purely quantum fluctuations are playing an increasingly 
important role and quantum character of thermal fluctuations 
should be taken into account. How to do it is not a simple task 


The stationary solution of this equation reads 

P(X) = N 0 D-' ix) expt-'P^x)], (25) 

where the generalized thermodynamic potential takes the form 

Aix)=pVix)-^[V'ix)] 2 , (26) 

We emphasize that the stationary distribution describes an equi¬ 
librium state, but it is not a Gibbs state. Remember that the 
Gibbs state is correct in the limit of a weak coupling of the sys¬ 
tem with thermostat. The Smoluchowski limit corresponds to 
the strong coupling regime. 

3. Discussion of results 

For the ring without a quantum dot, our model reproduces 
experimental data both for the diamagnetic and paramagnetic 
response in the vicinity of zero magnetic field E). For the 
ring with a quantum dot, we have not found experimental data. 
Therefore, our work could inspire experimentalists to design 
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Figure 2: The stationary probability distribution P(x ) of the dimensionless magnetic flux x in the classical Smoluchowski regime for the external magnetic flux 
x e = 0 and four values of the phase Sp of the transmission coefficient. In each panel there are four curves which correspond to different values of the amplitude 
of the transmission coefficient: tp = 0 (green dashed-dotted), 1/3 (black dotted), 2/3 (blue dashed) and 1 (red solid). Other parameters are: a = 0.1, ko = 1 and 
T 0 = 0.2. 



Figure 3: The stationary probability distribution P(x) in the quantum Smoluchowski regime for x e = 0 and four values of the phase Sp of the transmission coefficient. 
In each panel there are four curves which correspond to selected values of the amplitude of the transmission coefficient: tp = 0 (green dashed-dotted), 1/3 (black 
dotted), 2/3 (blue dashed) and 1 (red solid). The remaining parameters read a = 0.1, ko = l, To = 0.2, To = 0.001 and e = 100. 
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experiments and verify our theoretical predictions revealed be¬ 
low: the influence of the transmission coefficient t F and the 
phase S F of the quantum dot on stationary current-flux charac¬ 
teristics. 

The system has a 8-dimensional parameter space 
{x e , To,ko,a, To, e,t F , S F }. It would be difficult to carry 
out a comprehensive analysis and present current-flux char¬ 
acteristics for all possible sets of parameters. Therefore, for 
numerical calculations, values of the parameters a = 0.1, 
ko = l and Tq = 0.2 are kept fix. We include quantum 
corrections which are characterized by 2 parameters: To and 
6. Their physical meaning is explained in Refs. Il33ll34ll . The 
quantum dot is also characterized by 2 parameters: t F and 
S F and their impact is displayed below, they will be fixed at 
the value of To = 0.001 and e - 100. The similar analysis 
but for the pure metallic ring without the quantum dot is 
presented in our previous papers. The stationary solutions of 
the Fokker-Planck equation was addressed in Ref. [34] and the 
current-flux characteristics was investigated in Refs. [ 15 16]. 


3.1. Stationary states 

The stationary solution of the Fokker-Planck equations (]"3) 
and ( [23] ) is given by the steady-state probability distribution 
P(x) through the relations ( [18) and ( [25] ) without and with the 
quantum corrections, respectively. We consider the case = 0, 
i.e. when the external magnetic field is absent. In the case 
of classical thermal fluctuations, the Boltzmann distribution 
is depicted in Fig. [2] for four different values of the phase 
of the transmission coefficient S F = 0 (top-left), 1 (bottom- 
left), 7r/2 (top-right), n (bottom-right). For = 0, the prob¬ 
ability distribution is symmetric with respect to the reflection 
v —> -x. Moreover, it is invariant under the change of the 
phase P(x,S F ) = P(x,2n - S F ). Therefore below we con¬ 
sider the interval of the phase S F e [0,7r]. All four panels 
present the distributions for four different transmission coeffi¬ 
cient t F - 0,1/3,2/3,1. 

For full transmission (i.e. t F — 1) the distribution pos¬ 
sesses two local maxima for low valued phases, which reflects 
the bistability of the generalized thermodynamic potential 'Pc- 
This, in turn, means that in the steady-state the current can flow 
in two direction: clockwise or counterclockwise (but the aver¬ 
aged current is zero!). For the phase S F ^ 1 and full transmis¬ 
sion three local maxima can be found, with the most probable 
aside the local maximum around v = 0 (which denotes the zero 
current state). The additional local extrema, which doesn’t ap¬ 
pear in the S F —> 0 case, indicate the possible multi-stability. 
This means that again the self-sustaining persistent currents can 
appear without the applied magnetic flux and are more probable 
than the zero current state. For the moderate-to-high phases the 
local maximum of the probability distribution at v = 0 becomes 
the most protruding among all others located at more distant 
values of the flux x. It means that self-sustaining currents are 
difficult to induce. Moreover, the lifetimes of the induced cur¬ 
rents related to the remote from zero extrema are also expect to 
be relatively short El. 

As already stated, in the quantum Smoluchowski limit, the 
stationary solution ([25) describes the thermodynamic equilib¬ 


rium. It is not, however, the quantum Gibbs state, as we deal 
with the strong coupling to the environment. In this case, the 
probability distribution depends explicitly on the coupling of 
the ring with thermostat via the resistance R in the parameter 
To in Eq. The equilibrium stationary distribution with 

quantum corrections is depicted in Fig. [3] for the same set of 
the parameters as in the classical counterpart in Fig. [2] For 
the quantum corrections we set To = 0.001 and e - 100. This 
means that the difference between the quantum and classical 
fluctuations of the dimensionless magnetic flux is T = 0.0075. 
The corrected distribution display somehow magnified features 
seen in the corresponding classical cases: minima are deeper 
and maxima are more pronounced. 

3.2. Current-flux characteristics 

In previous papers icaiE). impact of quantumness of ther¬ 
mal fluctuations on the current-flux characteristics has been 
studied. In this section we will focus on influence of the quan¬ 
tum dot on such characteristics. For zero external magnetic 
flux, = 0, the averaged stationary current is zero. It fol¬ 
lows from the properties of the stationary distribution: it is an 
even function of x. The non-zero magnetic flux x e ^ 0 breaks 
the v-reversal symmetry and the non-zero averaged current can 
emerge. In Fig. [4] we depict the response of the metallic ring 
to the applied constant magnetic flux in the classical Smolu¬ 
chowski regime (i.e. for T = 0). It is worth to stress that the 
current characteristics for the amplitude t F — 1 and the phase 
S F = 0 of the transmission coefficient (top panel, red curve) 
represent the situation with maximal current. In other words, 
it is the same as the ring without quantum dot. The suppres¬ 
sion of the generated signal which comes with the reduction 
of the transmission amplitude seems to be the usual situation. 
For t F - 0 it is impossible to generate current in the ring. For 
the phase S F e (—tt/2, tt/ 2) the current response of the ring 
is paramagnetic for all non-zero amplitudes t F . In turn, for 
S F E (7i/ 2, 37r/2) the response is diamagnetic. Let us note the 
doubled period for the particular case S F = zr/2. The analysis 
for slightly lower or higher phases shows simple para- or dia¬ 
magnetic single-periodic structure of current-flux characteris¬ 
tics, respectively. 

We now address the issue of whether, and to which extent, the 
quantum nature of thermal fluctuations can influence transport 
properties. We thus show impact of quantum corrections on the 
current characteristics in Fig. [5] for the fixed quantumness pa¬ 
rameters To = 0.001 and e = 100. This figure is organized in 
exactly the same way as the previous one although the pecu¬ 
liarities are slightly different. For instance we cannot conclude 
here, that the maximal possible current amplitude is typically 
realized for t F - 1. In the classical Smoluchowski regime, the 
case t F - 1 is always the most optimal. With quantum correc¬ 
tions, it is intriguing to note that around S F = n/2 the maximal 
amplitude of the transmission coefficient does not provide max¬ 
imal current. In fact the current is weaker for t F — 1 than for 
t F — 2/3 or even when t F = 1/3, see Fig. 5. In fact one can ob¬ 
serve something similar to the transition from the paramagnetic 
to the diamagnetic state simply by changing the phase around 
S F = t:/2. This is displayed in Fig. [6] For the phases a little 
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tp —0 - tp — 2/3 






Figure 4: The stationary averaged current i in the classical Smoluchowski regime versus the external magnetic flux is presented for four different values of the 
phase Sf. In each panel there are four curves which correspond to selected values of the amplitude of the transmission coefficient: tp = 0 (green dashed-dotted), 
1/3 (black dotted), 2/3 (blue dashed) and 1 (red solid). Other system parameters read a = 0.1, ko = 1, Tq = 0.2. 


tp —0 ~ tp —2/3 






Figure 5: The stationary averaged current i in the quantum Smoluchowski regime versus the external magnetic flux is presented for four values of the phase Sp 
and four values of the amplitude of the transmission coefficient: tp - 0 (green dashed-dotted), 1/3 (black dotted), 2/3 (blue dashed) and 1 (red solid). Other system 
parameters read a = 0.1, ko = 1, Tq = 0.2, To = 0.001 and e = 100. 
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tp = 2/3 - tp = 1 — — t F =2/3 — — tp = 1 




Figure 6: The regimes of paramagnetic (/ > 0) and diamagnetic (i < 0) re¬ 
sponse in dependence of the phase Sp and the amplitude tp of the transmission 
coefficient through a quantum dot for an electron of the Fermi energy. The ex¬ 
ternal magnetic flux is fixed atx e = 1 /4. For the phase close to 7t/2 the current 
reversal is observed (upper panel). In bottom panel, the stationary averaged 
current in the quantum Smoluchowski regime is depicted as a function of the 
amplitude tp and in the vicinity of the phase 5p = tt/2. Four curves correspond 
to four values of the phase 8 p = 0.957r/2,7r/2, 1.017t/2, 1.05^/2. Other system 
parameters read a = 0.1, ko = 1, To = 0.2, To = 0.001 and e = 100. 

bit higher than tt/2, like one identify the classical picture - c.f. 
bottom panel on Fig. [5] 

As the next point of analysis we ask about domains of pa¬ 
rameters dp and t F where the current is paramagnetic and dia¬ 
magnetic, see Fig. [6] In the case of classical thermal fluctu¬ 
ations, the current is always of a paramagnetic type in the in¬ 
terval S F e (0,7t/ 2) and is always of a diamagnetic type for 
Sp e (tt/2, 3tt/3). In the case of quantum thermal fluctuations, it 
is not true: these intervals depend on the amplitude of the trans¬ 
mission coefficient. Nevertheless, the current is paramagnetic 
in a large interval around Sp = 0 and is diamagnetic in a large 
interval around Sp = n, and the transition point is in a small 
interval around 5p = n/ 2 . The type of response is more robust 
to changes in the amplitude of the transmission coefficient and 
more sensitive to changes of the phase around the value n/2. 

4. Summary 

This paper presents the influence of the quantum dot on trans¬ 
port properties of mesoscopic non-superconducting rings. The 
theory is constructed in the framework of the Langevin equation 
for magnetic flux dynamics. We have considered the case when 
the system is driven by classical thermal noise in the Smolu¬ 
chowski regime. The so named quantum Smoluchowski regime 


has also been studied. The stationary probability distribution 
both in ’classical’ and ’quantum’ case is depicted for zero exter¬ 
nal magnetic flux. The current-flux characteristics are analyzed 
in detail. The impact of parameters characterizing the quantum 
dot on the current has been addressed in this work. The phase of 
the transmission coefficient plays the crucial role in type of the 
response. In the ’classical’ case, its crossover value is fixed to 
Sp = 7t/ 2. Below this value, the current is paramagnetic while 
above this value the current is diamagnetic. For the ’quantum’ 
case, the response threshold depends on other parameters of the 
system, nevertheless it is located close to the value n/ 2. Finally, 
we would like to mention that recent progress in entirely novel 
experimental techniques makes the verification of our findings 
possible and we hope that our work will contribute to the de¬ 
velopment of effective control methods of transport properties 
in mesoscopic systems. 
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